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Quantum amplification is essential for various quantum technologies such as communication and 
weak-signal detection. However, its practical use is still limited due to inevitable device fragility 
that brings about distortion in the output signal or state. This paper presents a general theory that 
solves this critical issue. The key idea is simple and easy to implement: just a passive feedback of the 
amplifier’s auxiliary mode, which is usually thrown away. In fact, this scheme makes the controlled 
amplifier signihcantly robust, and furthermore it realizes the minimum-noise amplihcation even 
under realistic imperfections. Hence, the presented theory enables the quantum amplification to be 
implemented at a practical level. Also, a nondegenerate parametric amplifier subjected to a special 
detuning is proposed to show that, additionally, it has a broadband nature. 

PACS numbers: 42.65.Yj, 02.30.Yy, 03.65.Yz, 42.50.Lc 


I. INTRODUCTION 


The amplifier is clearly one of the most important com¬ 
ponents incorporated in almost all current technological 
devices. The basic function of an autonomous ampli¬ 
fier is simply to transform an input signal u to y = Gu 
with gain G > 1. However, such an amplifier is frag¬ 
ile in the sense that the device parameters change easily, 
and eventually distortion occurs in the output y. This 
was indeed a most serious issue which had prevented any 
practical use of amplifiers in, e.g., telecommunication. 
Fortunately this issue was finally resolved back in 1927 
by Black [T|, i; there are a huge number of textbooks 
and articles reviewing this revolutionary work, and here 
we refer to Refs. [1, lj|- The key idea is the use of feed¬ 
back shown in Fig. 1; that is, an autonomous amplifier 
called the “plant” is combined with a “controller” in such 
a way that a portion of the plant’s output is fed back to 
the plant through the controller. Then the output of the 
whole controlled system is given by 


y = g(^)u, 


G 

1 + GK' 


( 1 ) 


where K is the gain of the controller. Now, if the plant 
has a large gain G 1, it immediately follows that 
G(*^) Ki 1/K. Hence, the whole system works as an am¬ 
plifier, simply provided that the controller is a passive 
device (i.e., an attenuator) with K < 1. Importantly, a 
passive device such as a resistor is very robust, and its 
parameters contained in K almost do not change. This is 
the mechanism of robust amplification realized by feed¬ 
back control. Note, of course, that this feedback archi¬ 
tecture is the core of an operational amplifier (op-amp). 

Surely there is no doubt about the importance of quan¬ 
tum amplifiers. A pertinent quantum counterpart to the 
classical amplifier is the phase-preserving linear amplifier 
(in what follows, we simply call it the “amplifier”). 
In fact, this system has a crucial role in diverse quantum 
technologies such as communication, weak-signal detec¬ 
tion, and state processing 0 i i E [H B in, 0. 



Controller 


FIG. 1: Classical feedback-amplification scheme: G is the 
gain of an autonomous amplifier, and K is the gain of a passive 
controller. 


In particular, recent substantial prog ress in both theory 
and experiments [n [n [13,0, [ii iM m mui, HI 
has further advanced this field. An important fact is that, 
however, an amplifier must be an active system powered 
by external energy sources, implying that its parameters 
are fragile and can change easily. Because of this pa¬ 
rameter fluctuation, the amplified output signal or state 
suffers from distortion [23, IH, [131 ■ As a consequence, 
the practical applicability of the quantum amplification 
is still severely limited. That is, we are now facing the 
same problem we had 90 years ago. 

To make the discussion clear, let us here describe 
the general quantum-amplification process. Ideally the 
amplifier transforms a bosonic input mode bi to bi = 
51^1+172^21 where 62 is an auxiliary mode, and the coeffi¬ 
cients satisfy |5ip —I52P = 1 from [61, foj] = 1. Hence, the 
output 61 is an amplified mode of 61 with gain > 1. 
A typical example of an amplifier is the optical nonde¬ 
generate parametric amplifier (NDPA), in which case gi 
and g 2 are frequency dependent as shown later. How¬ 
ever, note again that the system parameters, especially 
the coupling strength of the pumped crystal, cannot be 
kept exactly constant, and eventually the amplified out¬ 
put mode 61 has to be distorted. 

Now the motivation is clear; we need a quantum ver¬ 
sion of the feedback-amplification method described in 
the first paragraph. The contribution of this paper is, in 
fact, to develop a general theory for quantum feedback 
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amplification that resolves the fragility issue of quan¬ 
tum amplifiers. The key idea is simple and easy to im¬ 
plement, i.e., feedback of the auxiliary output mode 62 
through a passive controller to the auxiliary input mode 
62 ■ Indeed, it is proven that the whole controlled system 
possesses a strong robustness property against parame¬ 
ter fluctuations, which thus enables quantum amplifiers 
to be implemented at a practical level. This type of con¬ 
trol scheme i^ in gen eral, called the coherent feedback 
[H, [2^ [sO, IHl meaning that an output field is 

fed back to an input field through another quantum sys¬ 
tem without involving any measurement process; hence, 
an excess classical noise is not introduced in the feedback 
loop. Now note that the auxiliary output &2 has some in¬ 
formation about bi due to their entanglement, though 62 
is usually thrown away in the scenario of quantum am¬ 
plification. Thus, we have an interpretation that the pre¬ 
sented scheme utilizes the signal-recycling technique 
for reducing the sensitivity, unlike the conventional use of 
it for enhancing the sensitivity of the gravitational-wave 
detector. 


In addition to the above-described main contribution, 
some important results are obtained. First, we see that 
the controlled system reaches the fundamental quantum 
noise limit even if some imperfections are present in 
the feedback loop. This means that precise fabrication of 
the feedback control is not necessary, which thus again 
emphasizes the feasibility of the presented scheme. Next, 
this paper proposes a type of NDPA subjected to a spe¬ 
cial detuning that circumvents the usual gain-bandwidth 
trade-off in the amplification process. A drawback of this 
modified amplifier is that, as will be shown, it is very sen¬ 
sitive to the parameter fluctuation. The presented the¬ 
ory has a distinct advantage in such a situation; that is, 
this issue can now be resolved by constructing a feedback 
loop. Therefore, as a concrete application of the theory, 
this paper proposes a robust, near-minimum-noise, and 
broadband amplifier. 


Finally, note that there are a variety of quantum ampli¬ 
fiers considered in the literature such as an optical back- 
action evasion amplifier [s^; however, the schematic pre¬ 
sented in this paper is essentially different from all those 
modifications in the following sense. While those mod- 
ihed amplifiers have their own purposes for improving 
the performance or achieving the goal in some specific 
subjects (e.g. back-action evasion), the feedback scheme 
is a device-independent and purpose-independent funda¬ 
mental architecture that must be incorporated in all am¬ 
plifiers. In fact, in the classical regime, the “operation” 
part of an op-amp has its own purpose (e.g., differenti¬ 
ation and integration), but any op-amp does not work 
without feedback. 


II. MODEL OF PHASE-PRESERVING LINEAR 
QUANTUM AMPLIFIER 

Let us begin with a specific model: the NDPA. This 
is an optical cavity system with two internal modes oi 
and 02. They are orthogonally polarized and obey the 
following Hamiltonian: 

H = wiojoi -I- uj2a\a2 + — 01026 ^*““*), 

with A G M the coupling strength between the modes, 
oji the resonant frequencies of Oi, and 2a;o the pump 
frequency. Also, in the above expression the rotating- 
wave approximation is taken under the assumption 2a;o ~ 
uii 002- The system couples with a signal input bi 
and an auxiliary (idler) input 62 with strength k. Then, 
in the rotating frame at frequency ujq, the dynamics of 
the NDPA is given by the following Langevin equations 

[Uii, 113: 


dt 


J Q-i “h XQj^ - 

- Vnbi, 

(2) 

da\ 

dt 


J 0,2 Xo\ - 

- y/^bl, 

(3) 


where Ai = ui—ujq and A2 = 002—ujq are detuning. Also, 
the output equations (boundary conditions) are given by 

bi = y/Hai -I- 61, b\ = + ^2- ( 4 ) 

Now the Laplace transformation of an observable Xt in 
the Heisenberg picture is defined by 

nOC 

x{s) := / e~^*xtdt, 

Jo 

where Re(s) > 0 . Then the Laplace transforms of &i, 
etc., are connected by the following linear equations: 

h{s) = gi{s)bi{s) +g2{s)bl{s), 

. . (s — ^-f fAi)(s-f ^ fA2) — A^ —kA 

9 i{s) = --, 32 s = 7^, 

D{s) D{s) 

D{s) = -b ^ -f jAi^ “ *^2) - A^. 

The stability analysis can be conducted in the Laplace 
domain; that is, for the amplifier to be stable, all roots 
of the characteristic equations of the transfer functions 
(i.e., poles) must lie in the left-hand complex plane. In 
the above case, particularly when Ai = A2 = 0, the 
characteristic equation is D{s) = s^-I-ks-I-k ^/4 —A^ = 0 ; 
hence, /c ^/4 — A^ > 0 must be satisfied to guarantee the 
stability of the NDPA. 

The quantum-amplification process is described in the 
Fourier domain s = iui with ui the frequency; that is, we 
consider the linear transformation at the steady state, 
bi{iuj) = gi{iuj)bi{iuj) g2{iuj)bl{iuj). Note that gi and 
32 satisfy \gi{iu})\'^ — \g2{iuj)\‘^ = I for all w. In partic¬ 
ular, when Ai = A2 = 0, the amplification gain at the 
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FIG. 2: Coherent feedback configuration for the autonomous 
amplifier G (plant). dz,d 4 „d^, and de are unwanted noisy 
input fields. 


resonant frequency w = 0 is given by 


l 5 i( 0 )| 


+ 4 A 2 
|k 2 _ 4 A 2 |’ 


and it takes a large number nearly at the threshold A ~ 
k /2 — 0 . Thus, hi is in fact an amplified mode of hi with 
gain |gi|. 

The above example can be generalized; any phase¬ 
preserving linear quantum amplifier is modeled as an 
open dynamical system with two inputs and two out¬ 
puts. Let us represent the input-output relation in the 
Laplace domain as follows: 


■ &i(s) ■ 

. ^ 2 ( 5 ) - 
G(s) 


Gii(s) Gi 2 (s)' 
G2i(s) G22(s)J’ 


( 5 ) 


where bi{s) is the Laplace transformation of bi, etc. The 
transfer function matrix G(s) at s = iuj (i.e. the scatter¬ 
ing matrix) satisfies 

\Giiiiuj)f - |Gi2(M|2 = IG 22 MI" - \G2i{iuj)f = 1, 
G 2 i(iw)G*]^(ia;) — G 22 (*w)G* 2 (fw) = 0 Vw. ( 6 ) 

Thus, |Gii(ja;)| represents the amplification gain. 


III. THE QUANTUM FEEDBACK 
AMPLIFICATION 

A. Feedback configuration 

Our control scheme is based on coherent feedback; that 
is, the controller is also given by a quantum system and 
is connected to the plant through the input and output 
fields. Note that, if a measurement process is involved 
in the feedback loop, it inevitably introduces additional 
noise. Now we take a passive system (e.g. a beam splitter 
and an optical cavity) as the controller, with two inputs 
&3,&4 and two outputs 63,^4; note that a single-input 


and single-output passive system has a gain equal to 1 
and thus does not work as an attenuator. We represent 
the input-output relation of this system in the Laplace 
domain as follows: 


llis) 

. blis) . 

K{s) 


bl{s 

. ^4(5) . 


Kiiis) Ki 2 {s) 
K2i{s) K22{s) 


( 7 ) 


Here, the creation operator representation is taken to 
make the notation simple. Because of the passivity prop¬ 
erty, the transfer function matrix K(s) is unitary in the 
Fourier domain; i.e., K{iu}yK{iuj) = I holds for all uj. 

We now consider connecting the controller to the plant. 
But unlike the classical case, where both the plant and 
the controller can be a single input-output system and 
arbitrary split or addition of signal is allowed, designing 
a feedback scheme in the quantum case is not trivial. For 
example, we could divide bi into two paths by a beam 
splitter and use one of them for feedback purpose, but 
in this case the resultant whole controlled system is not 
a minimum-noise amplifier. Instead, this paper proposes 
the following feedback connection as shown in Fig. 2 : 

h = bs, 62 = h, (8) 


which is, of course, equivalent to b^ = 63 and b^ = b\. 
Note that in Fig. 2 practical unwanted noises c?3,..., de 
are illustrated, but these modes are ignored for the mo¬ 
ment. From Eqs. ®, (0 and ([ 5 ]), the whole controlled 
system, with inputs bi, b\ and outputs &i, 63, has the fol¬ 
lowing input-output relation in the Laplace domain: 


'bi{s) ' 


. blis) _ 



G^^\s) 


^(fb) 

'-"12 

r^(fb) 

*-"22 


(5) 

(s) 


■ 61(s) 

. blis) 


where 

= [Gii - K2i{GiiG22 - Gi2G2i)]/(I - K 21 G 22 ), 
G^^^ = (Gi2K22)/(1 - K 21 G 22 ), 

G^^^ = {G 2 iKii)/{l-K 2 iG 22 ). 

G 22 ^ = [ 7^12 + G22{,KiiK22 — K 12 K 21 )]/{1 — K2iG22)- 


The matrix entries satisfy the condition corresponding 
to Eq. ([6]), i.e., |G^^^(ia;)p — \Gf^\iuj)\'^ = 1 Vw, etc. 
Finally, as remarked in Sec. II, for the whole controlled 
system to be stable, the controller should be carefully 
designed so that all poles of G-^^(s) must lie in the left- 
hand complex plane, as demonstrated in Sec. V. 


B. Robust amplification via feedback 

We now focus on the output of the controlled system 
in the Fourier domain, i.e., 

bi{iuj) = G^^\iuj)bi{iuj) + G^^\iuj)b\{iuj), 
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and the amplification gain (ioj) | especially when the 

original gain |Gii(za;)| is large. Note that G^^^ looks 
somewhat different from the classical counterpart ©; 
hence, it is not immediate to see if can be ap¬ 

proximated by a function of only the controller. Nonethe¬ 
less, the analogous result to the classical case indeed 
holds as shown below. 

For the proof we use Eq. ([ 6 ]) (below, we omit the vari¬ 
able iur). First, from |G2i||Gii| = IG22IIG12I together 
with the other two equations, we have |Gii| = IG22I and 
IG12I = IG21I. Also G11G22 - G12G21 = G 22 /G*, holds. 
Here, in the limit |Gii| —>■ 00, it follows that 

IG11G22 — G12G21I IG22I 1 „ 

\^\ ■ 

This implies that (G11G22 — Gi2G2i)/|Gii| converges to 
zero in this limit. As a consequence, we have 


|g 1 ^)| = 


Gii/IGiil - i^ 2 i(GiiG 22 - Gi 2 G 2 i)/|Gn| 


l/|Gii|-iF 2 iG 22 /|Gii| 
Gn/|Gii| 1 


-if 2 iG 22 /|Gn| 


\Ko 


Hence, in the frequency range where the plant has a large 
gain |Gii(*w)| » 1 , the whole controlled system amplifies 
the input bi(iw) with gain (ja;)| ~ l/|Ar2i(*w)| > 1 . 
Therefore we obtain the desirable quantum robust ampli¬ 
fication method via feedback; that is, thanks to the fact 
that the passive controller is much more robust compared 
to the original amplifier, even if Gn changes while main¬ 
taining a large value, the whole controlled system carries 
out robust amplification with stable gain l/|A'2i|. 


C. Feedback gain synthesis 


Here we conduct a quantitative analysis on the robust¬ 
ness property, which provides a guideline for synthesizing 
the feedback gain K. To see the idea clearly, let us again 
consider the classical case ©■ Let AG be the fluctua¬ 
tion that occurs in the plant G; then the fluctuation that 
occurs in the whole controlled system G*^^^ is calculated 
as 


AG(fb) 


G-k AG G 

l-k(G-kAG)A: “ 1 + GK 


AG 

(l + GiF) 2 ’ 


which as a result leads to 


AG(f^) 1 AG 
G(fb) “ l + GK ’ 


Hence, the gain sensitivity to the unwanted fluctuation 
can be reduced by the factor 1/|1 -I- GK\ by feedback. 
Equation (| 9 ]) suggests to us not to design G and K sep¬ 
arately; rather, what determines the performance of the 
controlled amplifier is the loop gain GK. Actually, while 
the controlled amplification gain, G^*^^ ~ l/AT, can be 


made bigger by taking a smaller value of AT, we should 
not design a too small K such that GK k. 0 ; in this case, 
Eq. ([ 9 ]) yields = AG/G, and thus there is 

no improvement in the sensitivity. The so-called Bode 
plot developed by Bode (e.g., see Ref. 0 ) is a powerful 
graphical method for synthesizing K as well as G, and 
it is now the standard tool for general feedback circuit 
design. 

Now let us try to establish the quantum version of the 
above discussion. Note in advance that a straightforward 
calculation, like Eq. cannot be carried out in the 
quantum case, but nonetheless a similar useful equation 
for determining the controller parameter K21 is shown. 
First, due to |Gii| = IG22I and G11G22 — G12G21 = 
G22/GI1, we find 

^(fb) _ |Giip — Ar 2 lG 22 _ G22 G22 — K21 

“ " GIi(l - K21G22) " ^ ■ 1 - K21G22 ’ 

which thus leads to 

|^(fb)| _ IG22 - K2 i\ 

' |l-if 2 lG 22 r 

The fluctuation of the controlled gain is given by 

AIG*-*^! = 1(^22 + AG22)* - K2i\ _ IG22 - AI21I 
|1 — Ar2l(G22 + AG22)| |1 — Ar2lG22| 


Then, from the general relation |a: -I- e| ~ |a:| -I- {xe* + 
x*e)/2\x\ with x,€ G C and |e| <C 1 , the normalized fluc¬ 
tuation of the amplification gain of the controlled system 
can be explicitly calculated as 


AIGi^^l 


|G' 


(fb), 

11 I 


1 -\K 2 i\^ 

\G* 22 -K 2 i\^ 



G22 K21 

1 — K21G22 



• ( 10 ) 


Next, noting that IG22I = |Gii| » 1 and thus |iF2i| ^ 
IG22I, we find from Eq. (ITOl) that 


iGiri 


1 

IG22P 



G 


* 

22 


1 — K21G22 



Thus, from the general relation |Re(a;?/)| < |a;||?/| with 
x,y G C, it follows that 


A|G^^^| 1 _ IAG22I . , 

~ |1 — Ai 2 lG 22 | IG22I 

Equation ©I) has a similar form to Eq. ([ 9 ]) and indeed 
provides us a guideline for feedback design. That is, as 
in the classical case, the balance of K21 and G22 deter¬ 
mines the ability of the controlled system to suppress the 
fluctuation (note that IAG22I > AIG22I = A|Gii|). In 
particular, we now deduce a similar conclusion as in the 
classical case; if we choose a too small K21 in addition to 
|Gii| » 1 such that the loop gain K21G22 is almost zero, 
then substituting the relation 


A|Gii| = AIG22I = Re(G 52 AG 22 )/|G 22 | 
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into Eq. (fTOl) we obtain = AjGn|/|Gii|. 

That is, in this case the fluctuation is not at all sup¬ 
pressed via feedback. To design an appropriate controller 
gain K21, the Bode plot of the loop gain K2i{ioj)G22{i^) 
is useful. 


IV. QUANTUM NOISE LIMIT 


Let us define the noise magnitude of b by 

(|A&| 2 ) := ^{AbAb'f + Ab'^Ab), Ab = b - (b). 


Then, through the ideal amplification process 61 = gibi + 
g2b\, the noise magnitude must be also amplified as 
(|A6ip) = |5ip(|A&ip) -I- |( 72 p/ 2 , where &2 is assumed 
to be in the vacuum. This implies the degradation of the 
signal-to-noise ratio: 


K&i)P . \{bi)\^ 

(|A&i|2) (|A5i|2)+yl (|A5i|2) 


Hence, the added noise 


(S/N). 


, |g2p |gip-i 

■ 2|gi|2 2|giP 

quantifies the fidelity of the amplification process [1, Q . 
In particular, in the large amplification limit \gi\ —>■ 00 
we find A —>■ 1/2, which is called the quantum noise limit. 

Up to now, the ideal setup is assumed, and the con¬ 
trolled system is driven by only the signal bi and the aux¬ 
iliary input 64, implying that it actually reaches the quan¬ 
tum noise limit in the large amplification limit. Hence, 
here we consider the following general case where some 
excess noise exists, as illustrated in Fig. 2 ; the plant is 
subjected to an unwanted noise da that enters into the 
system in the form bi = Gnbi + Gi2b\ + Giada; the 
controller is also affected by a noise ^4; furthermore, the 
feedback transmission lines are lossy, which is modeled 
by inserting fictitious beam splitters with additional in¬ 
puts ds and dg. Note that da,..., dg are all annihilation 
modes. Then the output of the whole controlled system 
has the form 


= G^^^&i+G'^^&^+G^^^da+G'^l’^dl+Glf d^+Gl^’dt 


Then, if the excess noises are all vacuum, the added noise 
in the feedback-controlled amplification process, denoted 
by , satisfies 


lim A^^^ = lim A^°'> 


1 I |Gi 3 p 

2 IGiiP’ 


( 12 ) 


where is the added noise of the plant. The proof of 
Eq. (na, including the detailed forms of G^^^, is given in 
Appendix A. This is a very useful result for the following 
reasons. First, in the large amplification limit the two 



Frequency [ K rad/s] 


FIG. 3: (a) Gain profile of the specially detuned NDPA 

without feedback, (b) Gain profile of the feedback-controlled 
system with parameters A = 5k and various p. 


added noises A^^'^ and are equal; as a consequence 
the second term in the right-hand side of Eq. m is a 
function of only the plant and cannot be further altered 
by feedback control. Hence, we have the following no- 
go theorem: If the original amplifier does not reach the 
quantum noise limit (i.e., lim|Cj^q_j,oo |Gi3|/|Gii| > 0), 
the controller can never remove this excess noise. On 
the other hand, notably, the imperfections contained in 
the controller and the feedback transmission lines do not 
appear in Eq. (HH). This means that a very accurate 
fabrication of the feedback controller is not necessarily 
required, which is a desirable fact from a practical view¬ 
point. Thus, if the original amplifier operates with the 
minimum added noise, the controlled system reaches the 
quantum noise limit as well even if some imperfections 
are present in the feedback loop. 


V. APPLICATION TO OPTICAL BROADBAND 
AMPLIFICATION 


In any practical situation, it is important to carefully 
engineer an amplifier so that it has a proper frequency 
bandwidth in which nearly constant amplification gain 
is realized. On the other hand, it is known in both 
the classical and quantum cases that, particularly for an 
amplifier with a single pole, the effective bandwidth be¬ 
comes smaller if the amplification gain is taken to be 
bigger. That is, there is a gain-bandwidth constraint. 
However, this constraint is not necessarily applied to a 
more complex amplifier with multiple poles. In fact, re¬ 
cently in Ref. , the authors propose a hybrid amplifier 
composed of two cavity modes and an additional opto¬ 
mechanical mode that circumvents the gain-bandwidth 
constraint. In this section, we study another system that 
is also free from this constraint, that yet does not need an 
additional degree of freedom. Then, the effectiveness of 
feedback is discussed, demonstrating its ability to make 
the system robust. 
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Plant 



FIG. 4: The NDPA (plant) and its coherent feedback; by 
simply feeding the auxiliary output 62 back to the auxiliary 
input 62 through just a beam splitter, the robustness property 
is drastically improved, as discussed in the main text. 

A. NDPA with special detuning 



Frequency [ K rad/s] 

FIG. 5: The upper blue lines represent the gain profile 

of the specially detuned NDPA without feedback, |Gii(iu;)|, 
while the lower red lines correspond to the controlled case 
with /3 = 0.1. In both cases, Ao = 5k. 


The plant system is the NDPA with dynamics ([ 2 ]), (|31) 
and output Q. Here we consider the ideal case where 
the unwanted noises .. ,dQ shown in Fig. 2 are not 
present. Without any invention, this system is subject 
to a gain-bandwidth constraint [^, but now let us take 
the specific detuning as Ai = A2 = A. The transfer 
function matrix of this system is then given by 


This expression shows that, as expected from the general 
theory discussed in Sec. Ill B, in the limit |Gii(za;)| —>■ 00 
(i.e., A —7> 00 ) the gain becomes |Gj^^(za;)| —>■ 1/(3 in a 
certain frequency bandwidth. To determine the attenua¬ 
tion level P, we need the stability condition; it is shown 
in Appendix B that, for all the poles of G^^^(s) to lie in 
the left-hand complex plane, the parameters must satisfy 


G(s) 


1 

(s -I- k/2)2 


— k^/4-I-zkA —kX 
—kX — k^/4 — inX 


The maximum gain is |Gii(0)| = -^/l -I- 16 A^/k^, which 
becomes larger by increasing A. Remarkably, this ampli¬ 
fication can be carried out without sacrificing the band¬ 
width. Figure 3 (a) shows the three cases corresponding 
to A = k,3k, 5k, all of which have the same effective 
bandwidth ~ k/10. A clear advantage of this system is 
in its implementability; that is, it is composed of only 
optical devices, and there is no need to prepare an auxil¬ 
iary system such as an opto-mechanical oscillator. Also, 
note that the system is always stable (the pole of the 
transfer function matrix is —k/ 2); that is, in a proper 
parameter regime such that the linearized model given 
by Eqs. @ and dS]) is valid there is no clear upper 
bound on A, in contrast to the standard NDPA which 
imposes |A| < k/2. 


B. The feedback effect 


This yields \I3\ < k/2|A| when /3^ <C 1; hence, let us here 
choose A = 5k, leading to |/3| < 0.1. Figure 3 (b) shows 
the gain |Gf;J’^(zw)| for the two cases /3 = 0.1 and /3 = 0.05 
together with the plot without feedback (i.e., (3 = 0). We 
then observe that the gain of the controlled amplifier be¬ 
comes smaller than that without feedback; in exchange 
for this reduced gain, the controlled amplifier obtains a 
great robustness property against the parameter fluctu¬ 
ation as demonstrated later. Note that a larger value of 
(3 (thus smaller amplification gain) induces a wider fre¬ 
quency bandwidth; hence, the controlled NDPA has the 
gain-bandwidth constraint. But the point here is rather 
that the gain and bandwidth can be easily tuned by just 
changing the reflectivity of the beam splitter. That is, an 
easily adjustable amplification can be realized, and this 
is also a clear advantage of feedback. 


Next, let us consider the feedback control of this am¬ 
plifier, again in the ideal setup. Here, as shown in Fig. 4, 
a beam splitter with transmissivity a and reflectivity (3 
is taken as a controller. This device has no internal dy¬ 
namics, and its transfer function matrix is constant: 


K{s) 


a (3 
(3 —a 


a, /3 S R. 


Thus, K 21 = (3 represents the attenuation level. The 
amplification gain of the whole controlled system is then 

(fb) _ (1 - ( 3 ) 8 ^ + (3ks - (1 -f / 3 )k ^/4 -k inX 
^ {1 — I 3 )s^ + Ks + (1 + ( 3 )k'^/ 4 : +i/SnX' 


C. Robustness property 

As repeatedly emphasized in this paper, the main 
strength of feedback is that the controlled system pos¬ 
sesses a robustness property. To see this, let us consider 
an imperfect case as follows. First, the device parame¬ 
ters are fragile; the coupling strength A fluctuates in such 
a way that A = (1 -I- 0.1eo)Ao, where Aq is the nominal 
value; similarly, the detunings Ai and A 2 can be slightly 
deviated from A, which is modeled by Ai = (1-1-0.001ei)A 
and A 2 = (1-1-0.001e2)A. Here (eo, ei, £ 2 ) are independent 
random variables subjected to the uniform distribution 
in [—1,1]. In addition to this fragility, we assume that 























7 








Y = 0.10K 

Y = 0.01K 



0.01 0.04 0.07 0.10 0.5 0.6 0.7 0.8 0.9 1.0 

Y/K ai 

FIG. 6: The added noise of the controlled amplifier with 
attenuation level /3 = 0.1 and that of the non-controlled one 
(i.e., /3 = 0) versus (a) 7 /k with fixed ai = 02 = 0.5, and (b) 
ai = 02 with fixed ')/k. In both figures, the red solid lines 
represent , while the blue dotted lines are at w = 0. 


the signal mode ai is subjected to optical loss, which is 
modeled by adding the extra term — 701/2 — to 

the right-hand side of Eq. with 7 the magnitude of 
the loss and ds the unwanted vacuum noise. The feed¬ 
back transmission lines are also lossy, which is modeled 
by Eq. (IMl) . 

The blue lines in Fig. 5 are 50 sample values of the 
autonomous gain |C?ii(ia;)| in the case Aq = 5k and 
ai = a 2 = 0.99. That is, in fact, due to the parameter 
fluctuation described above, the amplifier becomes fragile 
and the amplification gain significantly varies. Nonethe¬ 
less, this fluctuation can be suppressed by feedback; the 
red lines in Fig. 5 are 50 sample values of the controlled 
gain with attenuation level /3 = 0 . 1 , whose 

fluctuation is indeed much smaller than that of |Gii(ia;)| 
[ 5 ^ . (Note that, because the fluctuation of |Gj*|'^^(ia;)| 
is very small, the set of sample values looks like a thick 
line.) That is, the controlled system is certainly robust 
against the realistic fluctuation of the device parameters. 

D. Added noise 

Finally, let us investigate how much the excess noise 
is added to the output of the controlled or noncontrolled 
specially detuned NDPA. Again we set Aq = 5k, and 
the feedback control is conducted with attenuation level 
P = 0.1. Also, the same imperfections considered in the 
previous subsection are assumed; that is, the system suf¬ 
fers from the signal loss (represented with 7 ) and the 
probabilistic fluctuation of the parameters (A, Ai,A 2 ); 
furthermore, the feedback control is implemented with 
the lossy transmission lines (represented with ai and 02 ). 

With this setup Fig. 6 is obtained, where the red solid 
lines are sample values of the added noise at the center 
frequency w = 0 for the controlled system, given by 
Eq. (lAllI) . while the blue dotted lines represent those of 
the non-controlled system, given by Eq. (IA3I) . (In 
the figure it appears that six thick lines are plotted, but 
each is the set of 50 sample values.) Figure 6 (a) shows 

^(o) 

versus the signal loss rate 7 /K, where for 


the controlled system we fix a\ = = 0.5 (that is, 

the feedback transmission lines are very lossy). Also, 
Fig. 6 (b) shows the added noise as a function of a\ = 02 , 
with fixed signal loss t/k. 

The first crucial point is that, in both figures (a) and 
(b), A^^^ and are close to each other. This is the 
fact that can be expected from Eq. m, which states that 
_/^(o) coincide in the large amplification limit. It 
is also notable that, for all sample values, A^^'^ is smaller 
than A^°^ ; in other words, the feedback controller re¬ 

duces the added noise, although in the large amplification 
limit this effect becomes negligible as proven in Eq. m- 
Another important feature is that, as seen in Fig. 6 (a), 
the signal loss 7 is the dominant factor increasing the 
added noise, and the feedback loss 1 — ai(= 1 — 02 ) does 
not have a large impact on it, as seen in Fig. 6 (b). As 
consequence, when 7 is small, the controlled amplifier 
can perform amplification nearly at the quantum noise 
limit 1 / 2 , with almost no dependence on the feedback 
loss; this fact is also consistent with Eq. ( 121 ). 

In summary, the specially detuned NDPA with feed¬ 
back control functions as a robust, near-minimum-noise 
(if 7 <C k), and broadband amplifier. 

VI. CONCLUSION 

The presented feedback control theory resolves the crit¬ 
ical fragility issue in phase-preserving linear quantum 
amplifiers. The theory is general and thus applicable 
to many different physical setups, such as optics, opto- 
mechanics, superconducting circuits, and their hybridiza¬ 
tion. Moreover, the feedback scheme is simple and easy 
to implement, as demonstrated in Sec. V. Note also that 
the case of phase-conjugating amplification (^ . can 
be discussed in a similar way; see Appendix C. 

In a practical setting, the controller synthesis prob¬ 
lem becomes complicated, implying the need to develop a 
more sophisticated quantum feedback amplification the- 
ory, which indeed was established in the classical case 
@1 So, Sll • The combination of those classical approaches 
with the quantum control theory [H, IH, should ad¬ 
vance this research direction. Another interesting future 
work is to study genuine quantum-mechanical settings, 
e.g., probabilistic amplification [H, IH, [H, ■ Finally, 

note that feedback control is used in order to reach the 
quantum noise limit, in a different amplification scheme 
(the so-called op-amp mode) [SSi ; connection to these 
works is also to be investigated. 

This work was supported in part by JSPS Grant-in-Aid 
No. 15K06151. 


Appendix A: Proof of Eq. (12) 

First, we again describe the setup of the imperfect sys¬ 
tem depicted in Fig. 2, in a more detailed way. The plant 
system is subjected to an unwanted noise da, such that 












the input-output relationship is given by 
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The transfer function matrix in this case satisfies 


|Gll|^ — 10121 ^ + |Gl 3 p = |G' 22 |^ — |G 2 i|^ — |G 23 |^ = 1 , 
G2iGn - G22Gt2 + G 23 Gt 3 = 0 , (A 2 ) 


for all s = iu!. Thus, the added noise of the plant system 
is given by 
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This leads to 
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(A 4 ) 


where the second term is assumed to exist. Note that in 
a more general setup some creation noise modes (e.g. dg) 
can be contained in Eq. but this modification does 

not change the conclusion. The controller also contains 
an unwanted noise field in the following form: 
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(A 5 ) 

The point here is that only the creation modes appear 
in Eq. (IA 5 () . unlike Eq. dAH) that involves both creation 
and annihilation modes; this is indeed due to the pas¬ 
sivity property of the controller. Finally, the transmis¬ 
sion lines (optical fields) for feedback are assumed to be 
lossy. This setting is modeled by inserting two fictitious 
beam splitters; the beam splitter in the output field 62 
has transmissivity ai and reflectivity di, and also the 
beam splitter in the output field 64 has transmissivity 0:2 
and reflectivity S2 (we assume at, Si € R without loss of 
generality). Then, the coherent feedback connection is 
represented by the following relations: 


bl = aibl + Sidl, bl = a 2 bl + S 2 dl. (A6) 
Combining Eqs. (lATI) . (lAsI) . and (IA6I) . we end up with 


61 = G(‘^)6i+G(^)6i+G(^43+G(J4l+Glf d^+Gfe^’^d^, 


where 

G^^^ = [Gii - aia2if2i(GnG22 - Gi2G2i)]/G, 
gJ?^ = [Gi 3 - aia2if2i(Gi3G22 - Gi2G23)]/G, 
G^^) = a 2 Gi 2 X 22 /G, G^^^ = a 2 Gi 2 K 23 /G, 

G^^) = a 2 < 5 iGi 2 i^ 2 i/G, g 1 ^’ = < 52 Gi 2 /G, 


and G = 1 — aia2d^2iG22- Note that the transfer func¬ 
tions satisfy at s = *0; 
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Here we derive some preliminary results that are used 
later. First, Eq. (jA 2 l) leads to IG21G44 — G22G);2p = 
|G23p|Gi3p; together with the other two equations, we 
then have 
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Similarly, 
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(A 9 ) 


Furthermore, we can prove that IG22/G11I —)■ 0 never 
hold in the limit of |Gii| —>■ 00 as follows. If IG22/G11I —>■ 
0 , this leads to G21/G11 —^ 0 and G23/G11 —?> 0 , and fur¬ 
thermore, IG12/G11I —^ 1 and IG13I —^ 1 from Eqs. (IA 2 D 
and (|A8I) : then using Eq. (IA 9 I) we have G23 —>■ 0 and 
accordingly G2iG3^ — G22GI2 0, which leads to a con¬ 
tradiction. 

Now we are concerned with the amplification gain 
|Gf 4 l ill Ihs limit |Gii| —>■ 00. It is given by 


|Gf 4 l = 
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The second term is upper bounded, because from 
Eq. m\i 
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and also IG22/G11I does not converge to zero. Therefore, 
we need Ar2i(G22/Gii) —>■ 0 (hence, K21 —>■ 0) to have 
the condition IGg^^l —>■ 00. 

Finally, the added noise in the controlled system is 
computed as follows: 
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where Eq. (jA 7 | is used; also note that all the noise fields 
are now vacuum. The third term is given by 
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Then from Eq. 


we find 


X. 


GiiG'22 — G12G21 


21 
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11 


< \K 2 i\^ + 


K- 


21 


G^ 

Gii 


0 , 


in the limit -^121(^22/011) —>■ 0 and -^121 —t 0 . Also 
|-f^2i (G13G22 ~ G12G23)/Gii I —>■ 0 holds due to Eq. (IA 9 I) . 
As consequence, the added noise in the limit |gJ^^ | —>■ 00 
is given by 


lim 


4 (fb) 


1 I |Gl3p 

2 |Gii|2- 


Hence, together with Eq. (IA 4 I 1 . we obtain Eq. m- 
The point of this result is that, due to the strong con¬ 
straint on the noise input fields, which is represented by 
Eq. (TAtI) . the added noise does not explicitly contain the 
terms that stem from the creation input modes <^4,^5, 
and dg. This is because of the passivity property of the 
controller (IA 5 I) and the feedback transmission lines (IA6I) 
that are composed of only the creation modes. 


method [ 4 l| is a very useful graphical tool that can deal 
with such cases, although an exact stability condition is 
not available. Another way is a time-domain approach 
based on the so-called small-gain theorem that 

produces a sufficient condition for a feedback-controlled 
system to be stable; the quantum version of this method 
will be useful to test the stability of the controlled 
feedback amplifier. 



Appendix B: Proof of Eq. (13) 


FIG. 7: The poles represented by the cross points between 
two curves. 


The stability of the controlled amplifier is guaranteed 
if and only if all the poles of lie in the left-hand 

complex plane. In our case, those are given by the solu¬ 
tions of the following characteristic equation: 

(1 — / 3 )s^ -I- Ks -I- (1 -I- / 3 )k ^/4 -I- ifdKX = 0 . 

In the standard case where the coefficients of the char¬ 
acteristic equation are all real, the Routh-Hurwitz crite¬ 
rion can be used for the stability test, but now the above 
one contains an imaginary coefficient. Hence here we set 
s = X iy, x,y € M., transforming the above equation to 

K G 2 k^/ 3 ^ —( 3 kX 

2 ( 1 -/ 3 )J ^ 4 ( 1 -/ 3 ) 2 ’ 2{1 -/3)x + k 

The poles are given by the intersections of these curves 
in the complex plane; Fig. 7 shows the case for / 3 A > 0 . 
Hence, for the poles to be left in the complex plane, the 
parameters must satisfy 


Appendix C: Phase-conjugating case 


The Hermitian conjugate of the second element of 
Eq. ( 3 ) is given by 62 = G2ib\ G22b2- That is, the 
output 62 is the amplified signal of the conjugated input 
bl, with gain IG21I; this is called the phase-conjugating 
amplification. The feedback control in this case is almost 
the same as for the phase-preserving amplification. We 
consider the ideal feedback configuration shown in Fig. 2 
(i.e., the noise fields ds,..., de are ignored) and now fo¬ 
cus on the auxiliary output 63 = -f (G2^^)*64. 

Then the amplification gain is evaluated, in the large am¬ 
plification limit IG21I —)> 00, as 




Kii 

I/IG21I-A21G22/IG21I 


Kii 

-A 2 ie*« 


/ 1 

1^211 Y |if2lP 




Considering the other case (i.e. / 3 A < 0 ), we end up with 
the stability condition (HSl). 

Note that, as demonstrated above, in general the 
stability analysis becomes complicated for a complex- 
coefficient or higher-order transfer function. The Nyquist 


In the first line of the above equation, we have used 
IG22/G21P = 1-1- I/IG21P —1; also the last equality 
comes from the unitarity of K, i.e., I-Riip + |-f^2iP = 1 - 
Therefore, when the original amplification gain is large 
(IG21I » 1), the controlled system works as a phase- 
conjugating amplifier with gain y'l/|-A2iP — 1 > 1. As 
in the phase-preserving case, this controlled gain is ro¬ 
bust compared to the original one |G2i|. 
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